Charm and Bottom Quark Masses from Perturbative QCD 
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Using a new result for the first moment of the hadronic production cross section at order O(a^), 
and new data on the J/ip and ip' resonances for the charm quark, we determine the MS masses of 
the charm and bottom quarks to be m c (m c ) = 1.295 ±0.015 GeV and mb(rh~b) = 4.205 ±0.058 GeV. 
We assume that the continuum contribution to the sum rules is adequately described by pQCD. 
While we observe a large reduction of the perturbative error, the shifts induced by the theoretical 
input are very small. The main change in the central value of m c is related to the experimental data. 
On the other hand, the value of mt is not changed by our calculation to the assumed precision. 
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The strong coupling constant, a s , and quark masses 
are the only input parameters of the QCD lagrangian. 
Since at high energy masses become essentially negligi- 
ble, a s can be determined separately to a good precision. 
On the other hand, the determination of the masses is a 
much more cumbersome problem, where the details of the 
interaction can hardly be overcome. It is, therefore, not 
surprising that it has attracted a lot of attention. Par- 
ticularly interesting are the charm and bottom quarks, 
where results can be derived from lattice, see e.g. 1|, 
or from different types of sum rules, see e.g. 0,0, 4|. 
In this Letter, we shall follow the approach of [5j, which 
requires, as input on the theoretical side, moments of the 
hadronic cross section. The first moment at 0{or s ) con- 
stitutes our main new result 29], which we use, together 
with new resonance data for J / tp and ijj' , to derive the 
most up to date values of m c (m c ) and rn^rfib) in the MS 
scheme. 

The use of QCD sum rules for the determination of 
the charm quark mass has originally been proposed in a 
series of papers gathered in [lfj . The same approach has 
subsequently been applied to the bottom quark 0. a 
specific feature of these early studies has been the use 
of higher moments of the hadronic cross section. The 
advantages of the lower moments have been studied in 
"l|- and used at order 0(a 2 ) in and subsequently in 
jj |7j . The latter analysis was based on the calculation of 
three-loop moments from 0] and the data of the cross 
section scan around the charm threshold given in [l4j| . 
An important conclusion of Q is that the first moment 
is best suited for the analysis, since it has the weakest 
dependence on the nonperturbative effects and details of 
the threshold region, leading naturally to the smallest 
error. 

Starting from 0(al), the perturbative cross section 
receives contributions from diagrams with a massless 
threshold, Fig. Fortunately, gauge invariance restric- 
tions on the photon-three-gluon vertex allow such dia- 
grams only at 0((q 2 /m 2 ) 4 \og(q 2 /m 2 )) 0,0]. Since we 




FIG. 1: Singlet diagrams with a three-gluon cut. 

will use the first moment exclusively, we can safely ne- 
glect them in our analysis and restrict the calculation to 
diagrams where a single heavy quark line connects both 
electromagnetic current vertices. 

Let us now remind a few definitions. The hadronic 
ratio is given by 



R(s) 



cr(e + e — ► had.) 



(1) 



where er(e + e — > ) = Ana 2 /3s. It can be related 

to the current-current correlator through 

R{s) = 127rlmII( g 2 = s + ie), (2) 

where 

{-q 2 9^+q^quMq 2 )=i J dxe^ x (0\Tj^x)j v (0)\0). 

(3) 

For a given heavy quark of mass m, the contribution to 
the correlator can be expanded as 



lbn z 



n>0 



Am 2 (fi) 



(4) 



We perform all our calculations in the MS scheme, and 
expand the C n coefficients in the strong coupling as 



(5) 



The coefficients of the expansion in the strong coupling 
constant are known exactly up to two-loops and all mo- 
ments can be derived from a generating formula |l7i llSl ] . 
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The first eight moments of the three-loop result have been 
given in |13| . Our new result is the four-loop correction 

to the first moment, i.e. the coefficient C 1 . We have 
obtained it by a direct Taylor expansion to second order 
in q 2 of the photon self-energy, followed by a reduction 
of the resulting four-loop single scale tadpoles using an 



approach inspired by the Laporta algorithm [19j , see also 
[20j . We have kept a single power of the linear gauge pa- 
rameter in the Oth moment (value of n(0)) to check for 
gauge invariance, but neglected it in the 1st moment to 
spare unnecessary computational complexity. 

(3) 

The purely bosonic contribution to C ± is given by 
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where Cp and are Casimir operators of the funda- 
mental and adjoint representations respectively. In the 
case of the SU(N) group C F = {N 2 -l)/2N and C A = N. 
Furthermore, l m — log(m 2 (/x)//u 2 ), I2 — log 2, Q are Rie- 
mann zeta numbers, and di = Li^(l /2) with LU(x) the 
polylogarithm function. The numerical constants Ni are 
defined by the e expansion of the vacuum integrals given 
in Fig. |21 with the understanding that the integration 
measure per loop is e e7E /i7r 2 ~ e J d d k 



JVi = 


+2369.669517745751, 


JV 2 = 


+9090.208679977050, 


JV 3 = 


-764.0948373358558, 


JV 4 = 


-4647.352454831194, 


JV 5 = 


-1.808879546208335. 



(7) 



The above numbers have been obtained by changing 






... + N 1 e 2 + N 2 e i ... + N-ie + N 4 e 2 JV 5 e° 

FIG. 2: Definition of the numerical constants TV; from the e 
expansion of vacuum integrals. 



the normalization of the high precision results of |21j . 
In fact, we took the values of all the master integrals 
needed from there |3(ij| . Relations between the numbers 
in Eq. Q discovered in have not been used in our 
analysis. We checked, however, the values of all integrals 
up to six lines with the help of Mellin-Barnes integral 
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representations and the package \24 



7^(3) 



massless quark species is given by 



The single fermion loop contribution to C x with n; 
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where T F is the trace of the fundamental representation, The double fermion loop contribution to cf has been 
which we take to have the standard value T F = 1/2. Our published before in 25]. We found complete agreement 
result does not depend on the number of heavy quark with that result, and reproduce it here only for complete- 
species of mass m. It should be understood, however, ness 
that all terms without ni come from the heavy quark 
loops. 
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Substituting the color factors and the numerical values 
of the constants, our complete result can be cast in the 
following form 



1.87882 - 2.79472 n t + 0.0961014 nf 
+18.4050 l m - 5.26881 n t l m + 0.163146 nf l m 



(10) 



-0.660174 l 2 m + 0.474989 n t l 2 m + 0.0215821 nf lf n 
+0.656790 l 3 m - 0.256790 m 4 + 0.0197531 nf l 3 m . 

In order to determine the mass of a given quark from 



our perturbative result we shall use the moments of the 
hadronic cross section defined by 



M n = J ^R(s). 



(11) 



With the help of a dispersion relation, this can be trans- 
lated into 



12tt 2 / d . , 



q 2 =0' 



(12) 
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If one compares the above with Eq. 
finally given by 



10} , then the mass is 



m{jj) = 



1 /9 



2 \4 



77 X 1/(2") 



al, 



cxp 



(13) 



As explained at the beginning, we are only interested 
in the first moment. The experimental value, Eq. (|1 1|) . 
is made out of three constituents, the resonances, the 
threshold region and the continuum. We note that tiny 
changes in the central value of a s (Mz) [111 entering 
trough the perturbative result used for the continuum 
have no impact on the final moments. We observe the 
same with respect to parts of the five-loop result |27j . 
Therefore, apart from the resonance data for J/ip and 
ip' 0, we adopt the moments from Note that an al- 
ternative set of experimental moments that would lead 
to a larger error estimate can be found in |(| • After 
correction for the resonances, we have 



arm 

(-^1 P ) bottom 



= 0.2087 ±0.0042, (14) 
= 0.004456 ± 0.000121. (15) 



Similarly to [5j , we first determine the mass of the charm 
quark at 3 GeV, and then run it down to m c (rn, c ). For 
consistency reasons the running is performed with three- 
loop accuracy, i.e. to 0(al). We use the package [2|| 
and obtain 

m c (m c ) = 1.295 ± 0.009 Q ± 0.0003 M ± 0.012 oxp GeV. 

(16) 

The subscripts denote the various sources of the error 
(variation of a s within error bars, variation of fi in the 
range 3 ± 1 GeV and variation of the first moment). 
The combined error given in the abstract is obtained by 
adding the errors in squares. We note that the four-loop 
result gives a shift of -2 MeV whereas the new resonance 
data -10 MeV. Finally, we observe an almost threefold 
reduction of the perturbative error. 

Performing a similar analysis for the b-quark, where 
the calculation is first done at 10 GeV and then the mass 
is run down to m^mj), we obtain 

m b (m b ) = 4.205±0.010 cts ±0.002; I ±0.057 oxp . GeV. (17) 

The only change in the error determination with respect 
to the case of the c-quark is connected to the variation 
of [i in the range 10 ± 5 GeV. We note, that our result is 
exactly the same as that of the three-loop analysis in || , 
apart from a reduction of the perturbative error by al- 
most an order of magnitude. We would observe a shift of 
the central value, if the caculation were performed start- 
ing from low values of fj, close to 5 GeV. This is contained 
in our error estimate from the variation of ji. As a fi- 
nal remark, let us mention that there is some discussion 
on the resonance contribution of Y(4S') and T(55) 0. 
Adopting the values suggested by the latter paper, we 



observe a small shift of a few MeV in the central value, 
well within our error estimate. 
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